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The accelerated expansion of the universe is ascribed to the existence of dark energy. Black holes accretion 
of dark energy induces a mass change proportional to the energy density and pressure of the background dark 
energy fluid. The time scale during which the mass of black holes changes considerably is too long relative to the 
age of the universe, thus beyond detection possibilities. We propose to take advantage of the modified black hole 
masses for exploring the equation of state w[z] of dark energy, by investigating the evolution of supermassive black 
hole binaries on a dark energy background. Deriving the signatures of dark energy accretion on the evolution 
of binaries, we find that dark energy imprints on the emitted gravitational radiation and on the changes in the 
orbital radius of the binary can be within detection limits for certain supermassive black hole binaries. In this 
talk I describe how binaries can provide a useful tool in obtaining complementary information on the nature of 
dark energy, based on the work done with A. Kelleher. 



1. INTRODUCTION 

Dark energy drives the universe into an ac- 
celerated expansion. A decade after its discov- 
ery, the nature and origin of dark energy remain 
as elusive as ever. The detection of time varia- 
tions in the equation of state w[z] of dark energy, 
given by the ratio w[z] — p/p of its pressure p to 
the energy density p, is notoriously hard to con- 
strain. A precise measurement of w[z] will guide 
our theoretical exploration and shed light into the 
type of dark energy, thereby helping us under- 
stand our present universe and its future evolu- 
tion. For these reasons, a precision measurement 
of w[z] is of the highest priority in cosmology. 
Our current bounds are derived by a combina- 
tion of data from cosmic microwave background 
radiation (CMBR), Large Scale Structure (LSS) 
surveys, baryon acoustic oscillation (BAO) and 
Supernova IA (SnlA) data. Current data analy- 
sis constrains the dark energy equation of state to 
-0.14 < 1 + w < 0.12 0. While this represents 
progress, we would like to know if 1 + w is posi- 
tive, negative or 0, and whether or not w changes 



with time. Each case corresponds to quite differ- 
ent types of this mysterious energy and, it leads 
to dramatically different predictions for the future 
evolution of the universe. 

Dark energy evolves slowly with time for the 
cases when 1+w is positive or negative. However, 
since the bounds on its equation of state are close 
to a pure cosmological constant, w = — 1, then 
unlike matter dark energy does not cluster. For 
these reasons a universe filled with dark energy 
can be reasonably assumed to be a background of 
a perfect cosmic fluid. T.Jacobson and other au- 
thors |14|6j took dark energy to be a cosmic 'fluid' 
and showed that the mass of a black hole changes 
due to accretion of this background 'fluid'. The 
mass of the black hole increases or decreases, de- 
pending on the sign of 1+w. This effect would be 
hard to observe in a single black hole, as the time 
scale for this phenomena is quite slow relative to 
the age of the universe. For example the evolution 
time scale t for a solar mass black hole is about 
10 32 yrs. We proposed in [T] to use the evolution 
of supermassive black hole binaries instead of sin- 
gle black hole accretion to probe dark energy. The 
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reason for our proposed method relies on the fact 
that, instead of measurements of the evolution 
time scale of a single black hole, the 'footprints' 
of w[z] for supermassive black hole binaries can be 
observed and tracked down through the modifica- 
tions introduced by dark energy accretion in the 
the orbital separation of the binary and its emit- 
ted gravitational waves. For certain binaries dark 
energy modifications can be detectable by gravi- 
tational waves or X-ray and radio measurements. 
The information obtained by these modifications 
should help increase our bounds on w[z]. The use 
of binaries for probing dark energy is a different 
approach from our current methods of large scale 
experiments, since the observation utilizes local- 
ized systems, thereby avoiding noise inherited by 
propagation of signals through the vast structures 
of the universe. 

2. Dark Energy Accretion by Black Holes 

2.1. Evolution of a Single Black Hole on 
the Dark Energy Background: Review 

Let us assume that dark energy corresponds to 
a perfect fluid with an energy density p, pressure 
p and equation of state w[z] as a function of the 
redshift z, related by p — w[z]p. The case of a 
pure vacuum energy would have p — —p. As long 
as w[z] is close to -1 and varies slowly with time, 
the metric solution from Einstein's equations is 
taken to be approximately that of a De Sitter 
geometry, i.e. for w = — 1. 

A particularly interesting case is that of a 
Schwarzschild black hole in the background of the 
dark energy cosmic fluid. This case was stud- 
ied in [14131617] and the fluid accretion flowing 
through the black hole horizon was solved analyt- 
ically. The perfect fluid energy-momentum tensor 
with equation of state p = wp is assumed to be 



= p(l + w)u»u v + wpg^ 



(1) 



where g^ v is the inverse of the Schwarzschild 
metric, g^ v = diag(—(\ — 2A//r),(l — 
2M/r)~ 1 ,r 2 ,r 2 sin 2 9), and for the fluid 4- 
velocity, u^u^ = —1. Integrating the energy- 
momentum conservation equation and the pro- 
jection of the fluid 4-momcntum into it leads to 
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Following Babichev et al [7J, these expressions 
are manipulated into one for t 2 Tq. Then, in- 
tegrating the conservation law over the volume 
within the event horizon yields 



M = -iirr 2 ^ = AitAM 2 [ Poo + p( Poo )} 



(4) 



This expression can be integrated to give M 
as a function of time, neglecting the cosmological 
time evolution of p^: 
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The timescale for the accretion of dark energy 
Eq. [5] is given by the parameter t, and is r = 
j-. — . , .-^r 1 — ;— 7 — vtt ■ For a black hole of a mass 

(4ttAM(0) [poo +p(Poo )] ) 

m = am s which is a times larger than a solar mass 
m s , the evolution timescale is roughly r = 10 32 /a 
years. This is usually much longer than the age of 
the universe, thus beyond observational feasibil- 
ity. But the detection possibilities of the dark en- 
ergy accretion improve dramatically for the case 
of black hole binaries. The reason is that the 
evolution of the black hole binaries in the back- 
ground of dark energy is different from that of 
a single hole. Therefore, as described below, de- 
tection of the modifications to the orbital radius 
and the emission of gravitational radiation from 
these binaries in the background of the dark en- 
ergy fluid is within reach, [1] . 

2.2. Evolution of Black Hole Binaries in 
the Background of Dark Energy 

The Hulse- Taylor effect predicts a decrease of 
the orbital radius of the binary due the energy lost 
by the emission of gravitational radiation. This 
prediction has been succesfully tested. The evo- 
lution of black hole binaries in the background 
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of dark energy is modified relative to the Hulse- 
Taylor effect, due to the accretion of dark energy 
by the stars in the binary. Since the change in the 
mass of the stars, EqfSJ has a direct dependence 
on the parameters of the dark energy being ac- 
creted, specifically on w[z], then this information 
is carried out on the amplitude and the power 
of gravitational radiation produced PQ by the bi- 
nary. The dark energy accretion also imprints 
modifications in the orbital separation and merg- 
ing time. 

The modifications induced from the back- 
ground dark energy fluid onto the evolution of 
the binary, namely on the frequency uj of gravi- 
tational waves emitted by the binary and, on the 
orbital separation R, can be derived as follows: 
the change of the gravitational energy of the bi- 
nary is equal to the power lost due to gravita- 
tional radiation 5 : 



d ( lmim 2 \ ., 
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This is a temporal equation. Power losses 
via emission of gravitational radiation drives the 
binary's configuration to a new gravitational 
equlibrium separation. As a result the orbital ra- 
dius R decreases and eventually the stars inspiral 
and merge. For binaries immersed in the dark en- 
ergy fluid, it should be noticed that the point of 
gravitationally stable configurations of the binary 
at each moment is now driven by two effects: the 
usual loss of energy via gravitational waves emis- 
sion; and, the stars changing mass (leading to a 
change of the gravitational energy of the binary) 
due to dark energy accretion. The masses of the 
two black holes in the binary, mi and 7712, are in- 
creasing or decreasing with time, EqlSJ depending 
on (1 + w) being positive or negative. Therefore 
the two terms that induce changes in the orbital 
radius and period, namely: modifications due to 
dark energy accretion, and modifications due to 
energy losses from gravitational radiation com- 
pete with each other and determine the evolution 
of the orbit. 



Without loss of generality, this expression can 
be algebraically simplified by taking the two stars 
to be of equal mass, m,\ — m 2 = m with m = 
m(0). Replacing the constant mass of a black 
hole with the new expression, the mass rate of 
change of the black holes induced by the dark 
energy accretion from Eq|5j leads to a differential 
equation for the evolution of the binary R(t). 
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where the parameter L denotes L = 2G i mQT ■ This 
parameter contains all the modifications induced 
by the dark energy background and the new mod- 
ification terms due to dark energy to the orbital 
radius can be tracked down from all the terms 
in EqsElthat contain L. It should be noted that 
L ~ (l+w) is positive for quintessence type fluids; 
L is negative for phantom type fluids (l + w < 0); 
and, it becomes identically zero for w = — 1. In 
the latter case, all the modifications due to dark 
energy on the orbital radius vanish. The sign of 
L =/= determines the behaviour of the orbit, i.e. 
whether it grows or decreases with time for the 
cases when the 'L — terms' dominate over the 
conventional gravitational waves term. 

An approximate solution to Eql9] for the case 
when dark energy changes adiabatically is [I] 
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The fourth term in the expression Eq[9] corre- 
sponds to the conventional Hulse- Taylor term. 
As a consistency check, the Hulse- Taylor equa- 
tion is recovered in the limit when our universe 
approaches a DeSitter geometry w — > —I, (i.e. 
L — > 0). The Hulse- Taylor term describes the 
changes in the orbital radius that result for the 
energy losses of a binary from the emitted grav- 
itational radiation. The terms proportional to L 
are the new modification terms to the evolution of 
black hole binaries. They account for the effects 
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of dark energy accretion by the system. It should 
be noticed that one of the dark energy modifica- 
tion terms is of opposite sign to the Hulse- Taylor 
term; and, the type of dark energy with (1 + w) 
positive or negative leads to different types of evo- 
lution for the binary. An analysis of the solution 
for R[z] shows that the dark energy terms can 
dominate the evolution of the orbital radius, Eq|9j 
for certain cases of supermassive black hole bina- 
ries with large separation, quantified below. The 
different time evolution of the two terms, dark 
energy and the conventional Hulse- Taylor one, in 
Eq|9] on R and to, especially the linear depen- 
dence of R on the equation of state of dark energy 
1 + to, allow us to discriminate among the modi- 
fications to the orbit induced by dark energy and 
for probing the dark energy equation of state w[z] 
by observing the rate at which the orbit change 

it. 

In order to quantify the analysis of the above 
expression Eq[9]and discuss the interplay between 
the dark energy and Hulse- Taylor types of mod- 
ifications in the orbital radius we can parame- 
terize the binary as follows: let the initial mass 
of the star, (before modifications due to accre- 
tion), be too = am s where m s is a solar mass 
and a a parameter; and the orbital radius be (3 
times larger than the Schwarzchild radius of each 
star R = 2m (3^- Then, the ratio of the two 
correction terms to the binary's orbit R in Eql^l 
the Hulse- Taylor correction due to the emission 
of gravitational waves (GW) , and corrections due 
to dark energy (DE) accretion (terms containing 
L), is 

GW correction 10 45 ^ ^ 

DEcorrection (2/3) 5 (l + w)a 2 ~ 

We can use EqflO] to quantify the classes of bina- 
ries for which the dark energy correction terms 
dominate over the Hulse- Taylor term. Clearly 
for large enough separations (/3) or masses of the 
black holes (a ~ 10 18 ) the corrections due to dark 
energy can dominate over gravitational radiation. 

Observing this effect via gravitational waves 
experiments, we need both parameters of the bi- 
nary to be such that they favor the dark energy 
corrections over the Hulse- Taylor corrections to 
R, while at the same time being within observable 



ranges of frequency windows. The frequency and 
amplitude of gravitational radiation from these 
systems are given by 
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where r is the distance of the binary from the 
observer. EqfT^] constrains the second parame- 
ter (j not to be too large. For example LIGO is 
designed to detect radiation around 150 Hz opti- 
mally, and with an amplitude greater than h = 
1(T 23 Hz 1 / 2 .® LISA can detect much lower fre- 
quency radiation, down to ~ 10~ 5 Hz.]§\. Since 
for both GW experiments, the binary radius f3 
can not be too large, then our requirement of 
Eg 11 01 for using binaries to probe dark energy via 
their modifications on the orbital radius R, can be 
fulfilled by considering supermassive black holes, 
a ^> 1, (see pQ for specific examples). 

2.3. Observation Candidates? 

There have now been observations of black hole 
binary systems, and we can assess the effects of 
dark energy accretion on these systems. One ex- 
ample is galaxy 0402+379, observed in 2007 with 
VLBA. It has parameters a = 210 s , 2/3 = 10 6 , 
and r — 10 26 m. This case, while not observable 
through gravitational radiation, is one that is sen- 
sitive to the sign of 1 + w, and it's merging time is 
either greatly accelerated (1,000 years instead of 
60,000 years) by the effect of dark energy accre- 
tion when 1 + w > 0, or it is driven apart when 
1 + w < 0. 

Another example is Radio Galaxy OJ287, ob- 
served in 2008 with VLBA [3D]. This system is 
more complicated, since one of the black holes is 
more massive than the other. The parameters for 
this system are Rq ~ 10 20 to, or 2/3 ~ 10 7 and 
r ~ 3.5Mly ~ 10 22 to. While it is straightfor- 
ward to treat this system with the above equa- 
tions for toi ^ ma, the point is well illustrated 
by a simpler system of comparable parameters. 
We will use a ~ 10 9 . This gives / ~ 10~ 9 £fz and 
h ~ 10~20. Again, this system falls outside of 
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observable ranges for gravitational radiation, but 
the merging time is shortened by three orders of 
magnitude when 1 + w > 0, and the black holes 
are pulled apart for 1 + w < 0. 



of modified gravity and to discriminate those from 
the dark energy models. The evolution of binaries 
on the background of modified gravity is currently 
under investigation. 



3. Conclusions 
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